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ABSTRACT 

We use ideal axisymmetric relativistic magnetohydrodynamic simulations to cal- 
culate the spindown of a newly formed millisecond, B ~ 10 15 G, magnetar and its 
interaction with the surrounding stellar envelope during a core-collapse supernova 
(SN) explosion. The mass, angular momentum, and rotational energy lost by the neu- 
tron star are determined self-consistently given the thermal properties of the cooling 
neutron star's atmosphere and the wind's interaction with the surrounding star. The 
magnetar drives a relativistic magnetized wind into a cavity created by the outgoing 
SN shock. For high spindown powers (~ 10 51 — 10 52 ergs s _1 ), the magnetar wind is 
super-fast at almost all latitudes, while for lower spindown powers (~ 10 50 ergs s _1 ), 
the wind is sub- fast but still super- Alfvenic. In all cases, the rates at which the neutron 
star loses mass, angular momentum, and energy are very similar to the correspond- 
ing free wind values (< 30% differences), in spite of the causal contact between the 
neutron star and the stellar envelope. In addition, in all cases that we consider, the 
magnetar drives a collimated (~ 5 — 10°) relativistic jet out along the rotation axis of 
the star. Nearly all of the spindown power of the neutron star escapes via this polar 
jet, rather than being transferred to the more spherical SN explosion. The properties 
of this relativistic jet and its expected late-time evolution in the magnetar model are 
broadly consistent with observations of long duration gamma-ray bursts (GRBs) and 
their associated broad-lined Type Ic SN. 

Key words: Stars: neutron; stars: supernovae: general; gamma-rays: bursts; stars: 
winds, outflows; magnetic field; MHD 



1 INTRODUCTION 

Observations of long-duration gamma-ray bursts (LGRBs) 
have demonstrated that they are associated with core- 
collapse supernovae (SNe ) and the death o f massive star s 
(|Wooslev fe Blooml 120061 : iDella Valid 120061 : IZhand l2007h . 
Two leading candidates for the central engine powering 
LGR Bs are a newly formed, rapid l y rotating magnetar 
(e.g.. lUsovl Il99l iThompsori Il994l : " I Wheeler et all hood : 
iThompson et al. 2004) or a black hole with an accretion 
disk (e.g., [MacFadven fe Wooslevl Il999h . In both of these 
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cases, it is crucial to understand how relativistic material 
generated by the central engine - which is believed to give 
rise to the GRB at large distances from the compact object 
- escapes from deep within the host star. The most likely 
possibility, suggested by both afterglow observations ("jet 
breaks"; lRhoadslll999r ) and GRB energetics, is that a col- 
limated outflow punches through the stellar envelope, pro- 
viding a channel out of which the relativistic material can 
escape. 

In the collapsar model, collimated outflows from GRBs 
are accounted for by jets produced by a centrifugally sup- 
ported accretion flow onto a central black hole, based on the 
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results of numerical simulations fe.g. JProga et al. 2003) and 
by analogy with X-ray binary and AGN jets. 

In the magnetar model, however, the origin of such col- 
limated outflows is less clear becaus e relativistic outflows 
do n ot efficiently self-collimate (e.g., iLvubarskv &; Eichlerl 
l200lh . Numerical simulations of outflows from isolated mag- 
net ars have clearly shown that a highly collimated outflow 
can only be formed at very early times (within a few seconds 
after core bounce) when the outflow is only mildly relativis- 
tic, while relativistic outflows can only be produced at later 
times when the latitudinal distribution of the energy-flux 
approaches that of the force- free solution, with most of the 
energy fl ux confined in the equat orial region and not along 
the axis (Bu cciantini et al1l2006l ). These results appear to 
be independent of the specific magnetic field configuration 
on the surface of the star. 

Energy considerations demonstrate, however, that the 
surrounding stellar envelope provides an efficient confining 
medium even for a very energetic proto-magnetar wind. 
Thus, the GRB outflow might be strongly affected by the 
interaction with the progenitor star (thi s is true in both 
the c ollapsar and magnetar contexts). iKonigl &; Granotl 
2002 ) suggested, by analogy to pulsar wind nebulae 



Begelman feLilll992h (PWNe), that the interaction of the 



wind from the spinning-down magnetar with the surround- 
ing star could facilitate colli mation. 

In two previous papers (B ucciantini et al. 1 I2007L Paper 
I; iBucciantini et al] 120081 Paper II), we have investigated 
the dynamics of this interaction, using a variety of simpli- 
fying assumptions. In Paper I we used axisymmetric thin 
shell calculatio ns to demonstrate that collimation c an in- 
deed occur (see lUzdenskv fe MacFadvenll2006l . 120071 for re- 
lated ideas based on force-free rather than inertially loaded 
outflows). In Paper II we used relativistic magnetohydro- 
dynamic (MHD) simulations to investigate the interaction 
of a relativistic magnetized wind with the progenitor star, 
using wind s with properties ex pected for proto-magnetars 



based on | Metzger et all 20071) By analogy with PW Ne 



(bag 

([Komissarov &; Lvubarsky||2004l ; [Del Zanna et al.| [2004), we 
focused on the low-magnetization regime in which most (but 
not all) of the magnetic energy of the wind is assumed to 
be converted into bulk kinetic energy. We found that the 
magnetar wind shocks on the surrounding (exploding) stel- 
lar envelope, creating a bubble of relativistic plasma and 
magnetic field inside the star (a "magnetar wind nebula"; 
MWN). Just as in PWNe, if the toroidal magnetic field in the 
bubble is sufficiently strong, the bubble expands primarily 
in the polar direction due to the "tube of toothpaste" effect 
- the magnetic hoop stress in the equator creates a pressure 
gradient along the axis which is not balanced by any restor- 
ing force, thus driving the flow preferentially in the polar 
direction. The nebula itself is ultimately confined by the in- 
ertia of the SN ejecta, to which little energy is transferred, 
in contrast to the pres sure confinement in traditio nal mag- 
netic tower models (|Sherwin fe Lvnden- Bell 2007) or con- 
finement by a pressurized coco on inside the progenitor star 
(jUzdenskv fe MacFadvenll2006l ). 

The results of Paper I and II, despite being obtained 
with different approaches (thin-shell semi-analytic calcula- 
tions vs. full relativistic MHD simulations) and utilizing dif- 
ferent assumptions, demonstrate qualitatively similar over- 
all dynamics. In both cases, however, in order to reduce 



the computational time of the calculations, the properties 
of the outflow from the central source were specified using a 
fixed p rescription, based on previous results for fre e-flowing 
winds (|Bucciantini et al.ll2006l ; iMetzger et al.ll2007h , and us- 
ing super-fastmagnetosonic injection. One concern with this 
approach is that, if the termination shock between the wind 
and the surrounding star moves within the slow, Alfven, or 
fast critical surfaces, the stellar envelope is causally con- 
nected to the central star, which may modify how the star 
spins down and loses rotational energy and angular momen- 
tum (invalidating the boundary conditions used at the cen- 
tral source). In order to properly address these issues it is 
necessary to follow the dynamics from the proto-magnetar 

atm osphere out to large ra dii. 

iKomissarov &; Barkovl (|2007h carried out such a sim- 
ulation and found that there is indeed some causal con- 
tact between the stellar envelope and the central star. In 
spite of this difference relative to our assumptions in Pa- 
pers I and II, the collimation of the outflow by magnetic 
hoop stresses was simil a r to w hat we have previously found. 
IKomissarov &; Barkovl ([20071 V s simulations were only for 
~ 200 ms after core bounce, however, and thus only apply to 
the early non-relativist ic phase of magnetar spindown, not 
to the late-time relativistic, and potentially GRB-producing, 
phase. 

In this paper, we build on our previous work by carry- 
ing out time-dependent axisymmetric relativistic MHD sim- 
ulations of the development of magnetar outflows and their 
propagation through a surrounding star. We follow the dy- 
namics of the proto-magnetar wind into the late- time rela- 
tivistic phase (for ~ 10 sec), and from the injection radius 
located at the proto-neutron star surface to a distance of 
a few progenitor stellar radii. Given the computational re- 
quirements of this problem, which extends over several or- 
ders of magnitude in radius, and which involves timescales 
ranging from 10 ~ 7 sec in the proto-neutron star atmosphere 
to rsj 10 sec for jet propagation through the host star, we 
have not carried out a detailed parameter study, but have 
instead limited our investigation to a few fiducial cases. 

As in Paper II, we assume that an outgoing SN shock 
has already created a central evacuated cavity and that the 
host star is spherically symmetric. However, unlike in Paper 
II, we make no direct assumptions about the properties of 
the magnetar outflow. We only set the physical conditions at 
the proto-neutron star surface; specifically, the density, tem- 
perature and the radial component of the magnetic field are 
imposed there. A transonic flow self- consistently develops, 
with a speed increasing with radius. The wind magnetiza- 
tion a = Q 2 & 2 /M (where Q is the rotation rate, $ the open 
magnetic flux, and M the mass loss rate) changes in time as 
the mass loss rate drops due to the changing temperature 
and density in the neutron star atmosphere (which change 
as a result of the decreasing neutrino flux during the neu- 
tron star's Kelvin- Helmholz phase). Given that the outflow 
is not injected with fixed properties, but is instead allowed 
to self-consistently develop, any feedback due to the inter- 
action with the progenitor star is properly modeled. As a 
result (and unlike in Paper II), we do not assume that mag- 
netic energy is converted into kinetic energy at large radii in 
the wind. On the contrary, the outflows simulated here are 
always reasonably highly magnetized, from the early-time 
(< a few sec) non-relativistic wind in which the magnetic 
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energy is comparable to the kinetic energy at large radii, to 
the late- time magnetically dominated relativist ic outflow. It 
is, however, worth noting that our calculations may not cap- 
ture various instabilities (e.g., 3D ones) that could dissipate 
magnetic energy, converting it into thermal energy and ul- 
timately bulk kinetic energy; this limit of low a is explored 
in Paper II. 

The remainder of this paper is organized as follows. In 
J2]we describe our numerical methodology. We then describe 
the results of our calculations, including their implications 
for the collimation of magnetar outflows in LGRBs, the spin- 
down of the central neutron star, and nucleosynthesis of 
shock-heated stellar material (33]). Finally, in we sum- 
marize our results and their implications. 



2 NUMERICAL SETUP 

All of the simulations were performed using ECHO, a shock- 
cap turing central-scheme for g eneral relati vistic ideal MHD; 
see iDel Zanna &; Bucciantinil ((2002) and iDel Zanna et all 
( 2003|, 120071 ) for a detailed description of the equations and 
numerical algorithms. 

The interaction of the magnetar outflow with the sur- 
rounding SN progenitor is investigated by performing 2D 
axisymmetric simulations on a spherical grid. The domain in 
9 is the first quadrant from # = 0to# = 7r/2, with reflecting 
boundary conditions for vg, and the magnetic field com- 
ponents Be and at the polar axis to enforce axisymme- 
try, and similar boundary conditions in the equatorial plane. 
The grid is uniform in the 6 direction with 100 cells. Given 
that we are studying a wide range of spatial scales, from 
the proto-magnetar atmosphere at ~ 10 6 cm to the outer 
edge of the star at ~ 2 10 10 cm, we have selected a loga- 
rithmic spacing in radius with an inner boundary located at 
r m i n = 1.5 10 6 cm and an outer boundary at r ma x = 6.5 10 10 
cm, and 100 cells per decade in radius. Zeroth order extrapo- 
lation is assumed at the outer boundary. The code is second 
order in both space and time, with a monotonized central 
limiter, chosen in order to resolve the large density jump 
between the lighter relativistic plasma inside the MWN, as- 
sociated with the magnetar outflow, and the heavier stellar 
envelope (the density can increase by a factor of - 10 4 ). We 
use a Schwarzschild metric to account for the gravity of the 
central proto-neutron star, which must be included to prop- 
erly drive a transonic wind from the neutron star surface. 
We do not include the self-gravity of the progenitor star; 
the typical dynamical timescale for the progenitor is longer 
than the duration of our simulations and so the gravity of 
the progenitor can be neglected to first approximation. 

Our previous results (Papers I and II) imply that it 
takes ~ 5 — 10 s for the magnetar outflow to emerge from 
the progenitor surface and accelerate into the circumstellar 
medium. For this reason, we follow the evolution of the sys- 
tem for ~ 10 s. Note that most of the times quoted in this 
paper are given in seconds after core bounce. Our simulations 
begin 1 second after core bounce, which is approximatively 
the time it takes for the proto-neutron star to contract to 
its final radius c± 15 km and for the neutrino driv en wind 
to develop inside the SN ejecta ([Scheck et al.ll2006h . In or- 
der to evolve the dynamics in the neutron star atmosphere 
at 1.5 10 6 cm, the CLF time-step must be < 3 10~ 7 s. This 



highlights the fact that studying the evolution of the system 
for the desired duration of — 10 s is very time consuming. As 
a consequence, we have been forced to make several simpli- 
fications in our treatment of the microphysics. In particular, 
we use an ideal gas equation of state with an adiabatic index 
of 4/3, which is appropriate for both the relativistic magne- 
tar outflow and for the radiation pressure dominated shocks 
that result due to the interaction with the surrounding star; 
the validity of this assumption is checked in fc|3.4l 

The effects of neutrino heating and cooling, which 
are responsible for driving the wind, are primarily con- 
fined to a very small region around the prot o-neutron star 
([Thompson et al.1 l200ll ; iMetzger etafl l2007h . We approxi- 
mate the effect of neutrino heating by including an isother- 
mal layer at small radii near the inner boundary. At the 
inner radius we fix the rotation rate Q and the value of the 
radial magnetic field B r , and we assume perfect conductiv- 
ity so that Eq — E^ = 0. The temperature, density, and 
radial extent of the isothermal layer have been adjusted to 
reproduce reasonably accurately the mass loss rate obtained 
in more sophisticated ID calcul ations that include a detailed 
treatment of the microphysics ([Metzger et al.ll2007h , and to 
still guarantee numerical stability (decreasing the tempera- 
ture leads to numerical instabilities, especially at high mag- 
netization). We found that a reasonable choice was for an 
isothermal layer extending from the stellar surface to 3 10 6 
cm with a pressure p such that pc 2 /p c± 100 (i.e., a sound 
speed c s ~ 0.1c). The density p decreases in time oc t~ 2 ' 7 
in order to reproduce the decline in mass loss rate for a 
free wind that is caused by the decreasing neutrino luminos- 
ity of the neutron star in t he first 10 sec after core bounce 
([Burrows &; Lattimerl flQSG) ; the radial extent and the sound 
speed (set by p/ pc 2 ) in the isothermal layer are independent 
of time. Given that the wind is magnetocentrifugally acceler- 
ated, the value of the pressure does not affect its asymptotic 
dynamics (so long as the wind is cold). 

The evolution of the magnetization a, and the mass, en- 
ergy, and angular momentum loss rates for a free wind are 
shown in Fig. [3 (dashed line), discussed in ^3] A comparison 
between the mass-loss rates obtained using the isothermal 
layer approach described here with a I D model including a 
full treatment of the neutrino physics ([Metzger et al.l l2007) 
is shown in Table [1] and in Fig. The agreement is very 
good. Note also that the methodology adopted here pro- 
duce s the correct latitudina l variation in the mass-loss rate 
( see iBucciantini et al.l [2006). with magnetocentrifugal sup- 
port enhancing the mass-loss rate at the equator by a factor 
of ~ 100 relative to the pole for a neutron star with a ms ro- 
tation period. Although we have emphasized the comparison 
to free- wind calculations as a way of calibrating our simpli- 
fied model, it is important to stress that our boundary con- 
ditions only specify the density, temperature, and magnetic 
field in the neutron star's atmosphere; they do not specify 
the mass, angular momentum, or energy loss rates, which 
are determined self-consistently by the wind's dynamics. 

As in Papers I & II, we use the 35 M© model from 
IWooslev et ail ([20021 ) as our progenitor star. The outer sur- 
face of the progenitor is located at 2.5 10 10 cm. We assume 
that the density outside the star falls off as r -2 as expected 
for a wind. We have previously verified that for the range of 
radii we simulate, our results are independent of the outer 
density profile (Paper II). In order to account for the ef- 
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Table 1. Comparison of the free- wind mass loss rates obtained usi ng the isothermal lay er approach described in J2] with the values 
obtained by solving the full neutrino-heated MHD wind problem (Metzger et al. 20071), both using a ID monopole magnetic field. 
The cases with Period = oo correspond to non-rotating neutron stars, and are indicative of polar outflows in our 2D simulations. 
Time is after core bounce and mass loss rate M is in M@ 



Time (s) 


Period (ms) 


Magnetic field (G) 


M (Metzeer et al. 2007) 


M (this paper) 


1. 


oo 


10 15 


1.5 x 10~ 4 


1.6 x 10~ 4 


1. 


1. 


10 15 


1.2 x 10~ 2 


1.6 x 10- 2 


1. 


3. 


10 15 


3.0 x 10- 4 


4.0 x 10~ 4 


2. 


oo 


10 15 


3.0 x 10~ 5 


3.1 x 10~ 5 


2. 


1. 


10 15 


2.4 x 10~ 3 


2.9 x 10~ 3 


8. 


1. 


10 15 


9.3 x 10- 5 


9.0 x 10- 5 


1. 


1. 


10 16 


1.5 x KT 2 


1.9 x KT 2 


2. 


1. 


10 16 


2.6 x KT 3 


3.0 x KT 3 



feet of a SN shock propagating inside the progenitor, the 
region between 10 9 and 2 10 9 cm is given an initial out- 
ward velocity corresponding to a total kinetic energy 2 10 51 
erg, similar to that used in the one-dim ensional explosion 
calculations of IWooslev &; Weaver! ([l995). This corresponds 
to the SN shock moving at c± 1.5 10 4 km/s, 1 second af- 
ter core bounce. Hydrodynamic studies of neutrino driven 
winds after successful core-collapse SNe show that as the 
SN shock expands inside the star a cavity is left at small 
radii, inside which the wind blows ([Thompson et al.l 12005: 
IScheck et all [2006). To approximate this, our initial condi- 
tion includes a cavity inside the progenitor with a radius 
of 10 9 cm, which is roughly the size of the collapsing pro- 
genitor and the location of the SN shock 1 sec after core 
bounce. The cavity is initially filled with low density plasma 
in hydrostatic equilibrium. The density is chosen to be low 
enough that the proto-magnetar wind can initially expand 
freely into the cavity. The initial conditions do not, how- 
ever, determine whether the wind can expand freely at later 
times; this depends on the subsequent evolution, which is 
self-consistently calculated. 

The boundary conditions used here remove many of 
the assumptions about the wind structure that were used 
in our previous papers, and allow the outflow to develop 
freely and self-consistently according to the physics of the 
interaction between the wind and the progenitor star. The 
magnetic field is assumed to be monopolar and to extend 
from the surface of the proto-neutron star into the SN pro- 
genitor. This formally corresponds to assuming a large fos- 
sil magnetic field in the progenitor. We expect, however, 
that our results would be very simil ar if the field is gener- 
ated in situ by a convective dynamo ([Duncan &; Thompson! 
1992). For example, outside the light cylinder (~ 10 7 cm 
for a millisecond rotator), the structure and properties of 
the wind depend only on the amount of open magnetic flux, 
and not on its structure in the closed magnet osphere, so 
that the monopole solution is reasonable. Moreover, even if 
the poloidal magnetic field threads the entire star, at 10 9 cm 
(the initial size of the SN ejecta) it is dynamically negligible, 
and is unable to exert any significant torque on the central 
compact object. For the typical values of a in our simula- 
tions, the amount of open magnetic flux in the wind implies 
that the dipolar magnetic field is a factor ^> 2 larger tha n 
the monopole values quoted here ([Bucciantini et al1l2006h . 

To explore a few reasonable models for magnetar cen- 



tral engines of LGRBs, we consider three ca ses, motivated 
by previous work (e.g., iMetzger et al. I l2007h : Case A is a 
1 millisecond rotator with a B r — 10 15 C surface magnetic 
field; Case B is a 1 millisecond rotator with B r = 310 15 G; 
and Case C is a 3 millisecond rotator with B r = 10 15 G. In 
all cases the rotation rate is kept constant in time. After the 
fact, we can assess that this is a reasonable assumption for 
both cases A and C, in that only a modest fraction of the 
rotational energy is lost during the 10 seconds of our sim- 
ulation (see Fig. [71 discussed below). For case B, however, 
keeping the rotation rate constant is not fully self-consistent, 
because of the higher energy loss rate. Case B is nonetheless 
a useful guide to the dynamics in the high a, high spindown 
power limit. 

Relativistic MHD codes can have numerical difficulties 
at sufficiently large a and/or B 2 / p. In case C, we found that 
we were not able to lower the mass loss rate, and thus in- 
crease a, at late times, but had to artificially limit a to be 
^12; this is true even for the free wind. Despite the fact that 
the magnetization is lower than in cases A and B, case C was 
less numerically stable. We suspect that this is because the 
slower rotation rate in case C leads to a smaller mass loss- 
rate and thus lower densities. This affects the stability of the 
free wind, in which B 2 / p grows sufficiently large to cause nu- 
merical difficulties at large distances. In addition, in Case C 
there are more complex interactions between the magnetar 
wind and the supernova ejecta (see £I3.2[) : the compression 
of the magnetic field at small radii leads to the formation of 
strong current sheets, that the code fails to handle properly. 
This happens only at late times, and is due to the compres- 
sion of the magnetic field close to the Alfvenic surface at 
~ 10 7 cm, much smaller than the typical size of the MWN 
(see Fig. [6] discussed below). 



3 RESULTS 

We first describe the results of cases A and B, which are 
qualitatively similar. Case C has a significantly lower spin- 
down power, and has a somewhat different evolution; this 
will be discussed separately at the end. 
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Figure 1. Evolution of the radial velocity (in units of c) in a magnetized bubble inflated by a magnetar with B = 10 15 G and P = 1ms 
inside a 35 Mq star (Case A), about 2 seconds after core bounce (1 second after the start of the simulation). The time difference between 
the left and right panels is c± 0.15 sec. During this time, the termination shock collapses to small radii, inside the fast surface, near the 
pole, but it remains well outside the fast surface at the equator; the fast surface is at ~ 6 10 6 cm at the equator and ~ 1.2 10 8 cm at the 
pole. Once the termination shock becomes significantly curved near the pole, turbulent eddies are created, as can be identified by the 
negative radial velocity. 



3.1 Cases A &; B: High Spindown Power 

We first provide an overview of the evolution of the system 
as a function of time. The density in the initial cavity left 
by the expanding SN shock is small enough that the wind 
from the proto-magnetar can easily sweep through the cav- 
ity in a time ~ R C avit y /V W i n d ~ 0.1 sec. As soon as the wind 
impacts the high density eject a it is forced to slow down to 
a speed of order the SN shock velocity ~ 0.03c, in a strong 
termination shock (TS). This shock leads to the formation 
of a hot magnetized subsonic bubble (the MWN) whose evo- 
lution depends on the magnetization of the wind and on the 
spindown power of the proto-neutron star. 

In a ID monopolar geometry it is well-known that the 
TS can only exist for a long time at significant distances 
from the neutron star in the limit of a weakly magne tized 
wind (known as the a limit in PWNe; see lKennel Sz Coronitil 
1984). For a > Vmwn/c (where Vmwn is the expansion 
speed of the MWN, set initially by the SN shock velocity in 
our problem) the shock collapses towards the neutron star 
on a timescale of order a sound crossing time, due to the 
compression of the toroidal magnetic field. We have carried 
out ID relativist ic MHD simulations of magnetar spindown 
inside a star and have verified that for proto-magnetars with 
rotation rates and magnetic field strengths comparable to 
those considered here (see Table [IJ, the TS indeed collapses 
back down to the neutron star after ~ 0.03 sec. The result- 
ing causal contact between the central neutron star and the 
surrounding progenitor causes the neutron star to spin-down 
significantl y faster than is pre dicted by free wind calcula- 
tions (e.g.. iMetzger et al.ll2Q07h . 

The evolution is, however, significantly different in the 
multidimensional problem considered here. As the wind from 
the central neutron star inflates the MWN, two competing 
effects occur. On one hand, magnetic field gets progressively 
compressed inside the bubble, causing the TS to recede to 
smaller radii, as in the ID case. On the other hand, the high 



pressure inside the MWN pushes out against the denser SN 
ejecta. For a magnetized wind, the pressure along the rota- 
tion axis is significantly larger than at intermediate latitudes 
or the equator, for reasons that we discuss in detail in Papers 
I and II. This leads to a preferential expansion of the MWN 
in the polar direction and the formation of an elongated bub- 
ble. As a result, the a limit is less severe in 2D than in ID: 
magnetic field is allowed to flow from the equatorial region 
to the polar region in order to establish the magnetohydro- 
static pressure distribution (Beg elman &; Lilll992n . For large 
spindown power (Cases A & B), the expansion of the MWN 
in the polar region is fast enough to relieve the compression 
of the magnetic field by allowing an escape channel along 
the axis. As a consequence the TS moves towards the cen- 
tral neutron star at a slower speed than in ID. Once the 
jet emerges from the central part of the star, the outflow 
through the polar channel compensates for the compression 
of the magnetic field and the TS stabilizes. 

Figure [T] shows the radial velocity (in units of c) at three 
different times for Case A, all approximately 2 second after 
core bounce; the three times shown are within ~ 0.15 sec 
of each other. Initially, the TS is outside the fast surface at 
all latitudes and its shape is rou ghly spherical. As discu ssed 
in previous papers about PWNe (|Del Zanna et al.ll2004f ) the 
shape of the TS can be understood qualitatively in terms of 
pressure balance between the wind and the nebula. At early 
times the wind is moderately relativistic and is collimated 
in the polar direction; the ram pressure is thus higher along 
the axis than at the equator. The nebula also has a higher 
axial pressure due to t he compressed toroidal magnetic field 
(Beg elman &; Li| [l992). These two effects roughly balance 
and as a result the TS is roughly spherical. However, at 
later times the proto-magnetar wind becomes progressively 
more magnetically dominated, due to the decreasing neu- 
trino flux from the proto-neutron star, and the associated 
drop in mass loss rate. The spindown power thus becomes 
progressively more equatorial. At the same time, the com- 
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pression of the toroidal magnetic field in the MWN tends to 
increase the pressure anisotropy in the bubble. As a result 
the TS becomes oblate and a cusp forms at the pole, as can 
be seen at the later times in Figure [T] As soon as the polar 
cusp forms the inclination of the termination shock causes 
the formation of vorticity in the post shock region; this leads 
to the formation of large scale eddies along the axis. Because 
of the angle between the wind and the shock, the cusp also 
strongly reduces the effective ram pressure exerted by the 
wind. 

At later times, the TS along the axis moves to much 
smaller radii, inside the location of the fast surface (which 
along the axis it at a much large r distance from the n eutron 
star than on the equator; e.g. jBucciantini et al.ll2006h . Once 
the TS reaches the fast surface it becomes a weak Alfvenic 
discontinuity and a sub-fast outflow is established in the 
polar region which puts the nebula in causal contact with the 
central engine (Fig. [I]). However in both cases A and B the 
shock in the equatorial region never reaches the location of 
the fast surface, for two reasons: first, the TS shock tends to 
be at larger radii in the equator because the nebular pressure 
in the equator is lower than at the pole while for high a the 
wind ram pressure is higher at the equator; and, second, the 
fast surface is closer to the central engine in the equator (~ 
60 km at the equator and ~ 10 8 cm at the pole). The result 
is a mixed outflow in which the flow is subfast (but super- 
Alfvenic) in the polar region but superfast in the equatorial 
region. We will discuss the implications of this causal contact 
for the spindown of the neutron star in M3.3I 

In both cases A and B the expansion in the polar region 
is fast enough that even at relatively early times, the MWN 
has already created a jet-like feature along the rotation axis 
that is expanding at speeds that are a reasonable fraction of 
the speed of light. Figures [2] & [3] show the later time evolu- 
tion of the density, pressure, and velocity, for cases A and B, 
respectively. The polar jet expands much more rapidly than 
the SN shock and by ~ 3 — 5 sec the jet is outside the surface 
of the progenitor star. Despite some quantitative differences, 
the late-time evolution of the MWN is quite similar to that 
found in Paper II, in which we injected a high 7, low a wind. 

As soon as the jet starts expanding into the lower den- 
sity circumstellar region it accelerates to high speeds. For 
computational reasons, we cannot follow the acceleration of 
the jet for more than a few stellar radii, and thus we cannot 
completely assess the asymptotic acceleration of the out- 
flow. For the range of radii that we do study, we find that 
the Lorentz factor in the jet at late times is smaller than 
the value of a at the base of the proto-neutron star wind at 
the same time (the maximum achievable 7), but larger than 
the 7 c± a 1 / 3 acceleration of the ID monopole. Specifically, 
in case A we find a Lorentz factor 7^4 near the axis of 
the jet when a ~ 10 in the free wind at t ~ 10 sec; in case 
B, 7 ro 10 — 15 when a ~ 60 (see the top panel in Fig. 2}. 
The magnetization of the jet itself at large radii is also not 
negligible: in case A the jet has a magnetic field that is close 
to equipartition (B 2 / 'j 2 (pc 2 + 4P) /3 r — 1), while in case B it 
is somewhat above equipartition (B 2 /^ 2 (pc 2 + 4P)(3 r c± 5). 
This largely accounts for the fact that the acceleration is 
not 100% efficient. Indeed, the late-time magnetic + kinetic 
energy in the jet at large radii is comparable to the late-time 
value of the magnetic energy in the magnetar wind at small 
radii. This is a non-trivial result since the MWN tends to 



accumulate the flow so that the highly magnetized late-time 
outflow could have been partially mixed with the earlier less 
magnetized outflow. 

Although the magnetar wind efficiently and rapidly cre- 
ates a collimated polar jet that allows relativist ic material 
to escape the host star, the overall interaction with the rel- 
atively spherical SN ejecta is much weaker; in particular, 
little of the magnetar spin-down power is transferred to the 
spherical SN ejecta. This is true even at late times when 
the energy flux in the free magnetar wind at small radii is 
largely equatorial. The reason is that the wind undergoes a 
termination shock and then escapes via the polar channel. 
The low equatorial pressure in the MWN leads to very little 
energy transfer to the ejecta near the equator. In case A, 
we do find that the MWN is able to partially compress the 
shocked SN ejecta in the first few seconds, but as soon as the 
jet develops the internal energy of the MWN escapes via the 
polar channel and the SN ejecta start to recede inside and 
partially recompress the MWN. However, in neither case A 
nor case B do we find any significant changes in the global 
properties of the SN shock inside the progenitor star due to 
the magnetar wind. 

To quantify this, Figure |4] shows, for cases A and B, 
the relation between the instantaneous power in the mag- 
netar wind at small radii, the energy flux in the relativistic 
core of the jet at large radii (defined to be the energy flux 
within 5° of the pole) and the energy flux in the wider angle 
wind at large radii (defined to be the energy flux within 20° 
of the pole) . Figure |4] demonstrates that essentially all of 
the energy injected by the central engine is carried away in 
the collimated jet with a significant fraction of the energy 
confined in the central relativistic core. Note that there are 
~ 10 — 20% fluctuations in the energy flux due to turbulence 
created in the curved TS near the axis; there are significant 
fluctuations in the Lorentz factor as well (top panel). It is 
also interesting to note that the energy flux in the jet and 
wind at large radii is on average a little higher than the en- 
ergy flux in the wind at small radii at the same time. This is 
because energy injected by the wind at earlier times can be 
stored in the MWN, and released later on. The MWN thus 
introduces a small delay between conditions in the wind at 
injection and conditions in the jet emerging from the star. 

The vorticity created at the curved TS in the polar 
region leads to the development of a layer where sausage 
modes grow unstable. It is reasonable to expect that in a 
more realistic (e.g., 3D) calculation these instabilities might 
lead to the dissipation of some of the toroidal field. In 2] we 
will discuss how these instabilities might modify the evolu- 
tion of the system. 

3.2 Case C: Lower Spindown Power 

We now discuss case C, which corresponds to a magnetar 
with P = 3 ms and B = 10 15 G (see Fig. [5]). This case 
differs from the previous two in two important ways: first, 
the typical spindown power is a factor of ~ 10 (~ 100) 
lower than in case A (Case B); and, second, the location 
of the light cylinder and thus the fast surface is at a larger 
distance. Because of the lower spindown power, the evolution 
of the MWN proceeds more slowly than in Cases A and B. 
Ten seconds after the launching of the wind, the MWN has 
barely emerged from the progenitor star. 
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Figure 2. Evolution of the magnetized bubble inflated by a magnetar with B = 10 15 G and P = 1ms inside a 35 Mq star (Case A). From 
left to right, the panels show: logio [density (g cm -3 )], logio [pressure (erg cm -3 c 2 )] and velocity (in units of c). From top to bottom, 
the snapshots are 2,4, & 7 seconds after core bounce. The radius of the progenitor star is 2.5 10 10 cm. By t ~ 5 sec (middle panel) the 
jet has escaped the progenitor star. 
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Figure 3. Evolution of the magnetized bubble inflated by a magnetar with B = 3 x 10 15 G and P = 1ms inside a 35 Mq star (Case 
B). From left to right, the panels show: log io [density (g cm -3 )], logio [pressure (erg cm -3 c 2 )] and velocity (in units of c). From top to 
bottom, the snapshots are 1.6,3, & 5 seconds after core bounce. The radius of the progenitor star is 2.5 10 10 cm. By t ~ 3 sec (middle 
panel) the jet has escaped the progenitor star. 
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Figure 4. Outflow and spindown properties for Cases A (left) and B (right). Upper panels: Lorentz factor in the jet along the axis as a 
function of radius, at the end of the simulation (t = 11 s after core bounce). Lower panels: Comparison of the energy flux in the magnetar 
wind at small radii near the neutron star (solid line), energy flux in the jet at large radii (dotted) and energy flux in the cocoon at large 
radii (dashed). The jet is defined to be within 5° of the axis, while the cocoon is within 20° from the axis. 



The evolution of the MWN is also qualitatively different 
from the higher spindown power cases. The lower pressure 
in the MWN implies that the nebula expands at a signifi- 
cantly slower speed. As a result, even though the pressure 
along the rotation axis is larger than that at the equator, 
the expansion of the polar jet is not fast enough to overcome 
the compression of the toroidal field in the nebula, and so 
the termination shock collapses to smaller radii (unlike in 
the higher spindown cases where the rapid polar expansion 
allows the system to avoid this fate). A sub-fast flow is first 
established in the polar region, but after about 1 second the 
termination shock contracts within the fast surface along the 
equator as well, resulting in a fully sub-fast outflow. The dy- 
namics in this case is qualita tively similar to that found by 
Komissarov & Barkovl {2007) although the absolute radial 
scales are larger. In £13.31 we will discuss the implications of 
the sub-fast outflow for the spindown of the proto-magnetar. 
As in the previous cases, vorticity is created downstream of 
the termination shock that generates turbulence which per- 
sists to late times. This turbulence gives rise to significant 
interpenetration of the wind and the SN ejecta, in the form 
of fingers of SN ejecta that are dragged down to small radii, 
as is shown in the zoom- in to small radii in Figure [6] At late 
times, the denser SN ejecta even compress the magnetar 
wind at radii close to the Alfvenic surface. 

We have attributed the differences between cases A/B 
and case C - namely a largely super-fast vs. fully sub-fast 
outflow - to the differences in spindown power and rotation 
rate. One concern might be whether these properties of the 
outflow are sensitive to the size of the "SN" cavity we initial- 



ize. To assess this, we have repeated case A using a cavity 
with an initial radius of 5 10 8 cm, a factor of 2 smaller than in 
our fiducial case shown in Figure [2] At very early times < 1 
sec the smaller cavity does modify the solution, leading to 
a fully sub-fast outflow because the termination shock is at 
smaller radii. However, as the system evolves and the MWN 
expands into the outgoing SN ejecta a super-fast outflow is 
again established in the equatorial region. After ~ 1 sec, the 
system relaxes to a configuration similar to that which we 
found starting with a larger initial cavity. This demonstrates 
that the primary physics determining whether the magnetar 
wind is super-fast or sub-fast is indeed the spindown power 
of the neutron star. 



3.3 Proto-Neutron Star Spindown Rate 

In all of our simulations, there is at least some fraction of the 
proto-neutron star outflow that is in causal contact with the 
surrounding stellar envelope. In the higher spindown power 
simulations (cases A and B), the polar region is sub-fast, 
while in the lower spindown power simulation (case C), the 
entire outflow is sub- fast. In all three cases the outflow is 
always super- Alfvenic. Because of this causal contact, it is 
possible that the spindown of the neutron star could be mod- 
ified from that due to a free wind. Understanding the torque 
on the central engine is important because it strongly influ- 
ences whether or not conditions suitable for producing GRBs 
can be achieved. For example, an increased torque due to a 
sub-fast outflow might lead to sufficiently rapid spindown 
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Figure 5. Evolution of the magnetized bubble inflated by a magnetar with B = 10 15 G and P = 3ms inside a 35 Mq star (Case C). From 
left to right, the panels show: logio [density (g cm -3 )], logio [pressure (erg cm -3 c 2 )] and velocity (in units of c). From top to bottom, 
the snapshots are 2, 6, & 11 seconds after core bounce. Distances are in units of 10 9 cm; the radius of the progenitor star is 2.5 10 10 cm. 
By t ~ 10 sec (bottom panel) the jet has just started to escape the progenitor star. 
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Figure 6. Density (Left panel) and poloidal velocity (in units of c; Right panel) in the inner 10 8 cm near the end of our simulation of 
Case C (B = 10 15 G and P = 3ms). The vertical solid line near a cylindrical radius of 170 km in the velocity plot is the location of the 
Alfven surface. Note that the dense and slow moving stellar ejecta has compressed and redirected the magnetar's wind almost down to 
the Alfven surface. This leads to the enhanced torque and larger spindown variability seen in Fig. [7] 
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Figure 7. Magnetization a and mass, angular momentum, and energy losses for cases A (left), B (center), and C (right). The dashed 
line represents the same quantities computed for the free wind; for case C, the free wind was numerically unstable at late times and so is 
not plotted after ~ 6 s - For case A, we also show the mass loss rate at 3 times (stars) from the calculations of iMetzger et al.l (|2007l ). 
who included realistic neutrino microphysics. 
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that most of the rotational energy has been depleted before 
the outflow becomes relativistic. 

Figure [71 shows the magnetization, mass loss rate, en- 
ergy loss rate and angular momentum loss rate for the cen- 
tral neutron star in our simulations (solid) compared to the 
corresponding results for free winds (dashed lines). Figure [3 
shows that for cases A and B, the energy and angular mo- 
mentum losses in the wind are essentially identical to those 
of a free wind. At first, this seems to contradict the fact that 
a sub- fast wind has a larger torque than a super-fast wind 
in ID calculations. However, for cases A and B the sub-fast 
outflow is present only in the polar region, while at high a 
most of the torque is exerted by the equatorial outflow. The 
polar outflow thus makes only a small contribution to the 
total torque even if it causally connected with the envelope. 

The close similarity between the torque in our simula- 
tions and that exerted by a free wind can also be under- 
stood in terms of the location of the critical surfaces. In ID 
there are 2 free parameters at the inner boundary (i.e. 
and v r ), and 2 corresponding critical points, the slow and 
fast magnetosonic points; the equations can be renormal- 
ized to eliminate the singularity at the Alfven point, which 
is automatically crossed once the solution passes smoothly 
through the slow and fast points. These two critical points 
correspond to two eigenvalues of the system: the mass loss 
rate, associated with the slow point, and the terminal Lorenz 
factor, associated with the fast point (the torque, which is 
also an eigenvalue, is determined by the Alfven point, and, 
in ID, is fixed once the other two are known). In 2D, how- 
ever, there are three free parameters at the inner boundary 
(Be, in addition to v r and B^) and 3 corresponding critical 
surfaces and 3 independent eigenvalues: the slow, fast, and 
Alfven surfaces. In order to modify the torque significantly, 
the location and shape of the Alfven surface must be modi- 
fied. For fast rotators like those considered here, this implies 
that the outflow must be modified inside the light cylinder. 
We find, however, that even if the termination shock moves 
inside the fast surface, the location of the Alfven surface 
remains relatively unchanged. 

The results for case C - in which the outflow is fully 
sub- fast - are somewhat different than those for Cases A/B. 
Although the torque and energy losses are nearly identical to 
that of the free wind at early times, at late times the energy 
losses tend to be more variable and can be ~ 30% higher 
than for the free wind (Fig. [7}. These large fluctuations cor- 
respond to when the post termination shock turbulence is 
able to compress the proto-magnetar wind down to radii 
that are close to the light cylinder (Fig. [6}. This level of 
fluctuations, while interesting from the point of view of the 
underlying physics of the problem, is unlikely to change the 
overall evolution of the system in a significant way. 



3.4 Post-Shock Nucleosynthesis 

Our calculations have assumed an ideal gas equation of state 
with an adiabatic index of 4/3, which is appropriate for a 
radiation dominated flow or for relativistic conditions. We 
have checked a posteriori that this assumed adiabatic index 
is consistent with the thermal properties of the MWN and 
jet. Figure [8] shows the temperature derived from inverting 
the relation: 
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some nucleosynthesis of Carbon and Oxygen will occur in mate- 
rial shock heated by the polar jet. 
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between the pressure and density in our simulations, at t = 
1.6 sec after core bounce. At almost all locations, the tem- 
peratures are high enough that radiation pressure dominates 
and thus our equation of state is self-consistent. 

The temperature of the gas is also important for un- 
derstanding nucleosynthesis during the SN explosion. For 
example, explosive n ucleosynthesis of 56 Ni requires temper- 
atures of > 510 9 K ([Wooslev et alJl20Q2h . Figure [8] shows 
that this temperature is not attained even at relatively early 
times in our simulations. The reason for this is the low den- 
sity of the outer stellar envelope. By the time the jet-plume 
emerges outside the SN shock the density of the progenitor is 
~ 10 4-5 g cm -3 . At these densities Ni production requires 
a shock moving at nearly the speed of light, significantly 
faster than what we find even for our most energetic explo- 
sion (case B). We do find, however, that ~ 10 _2 M© of high 
speed (v c± 0.1 — 0.2 c) Ne and O can be created, because 
these have lower threshold temperatures for successful ex- 
plosive nucleosynthesis. The synthesized mass is relatively 
low because only a small solid angle near the pole is shock- 
heated to sufficiently high temperatures by the jet. In the 
context of observed GRBs, our viewing angle is relatively 
on-axis, in which case the high velocity nucleosynthesis may 
be observable; high velocity O and Ne can also be produced 
by the jet blowing out stellar material that had be en pro- 
cessed during stellar evolution (Ma zzali et al . 2006). Even 
if the jet nucleosynthesis only contributes marginally to the 
total nucleosynthesis during the explosion, it might lead to 
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unique observable signatur es in the eject a at la te times (as 
may be the case for Cas A; Wheeler et al.ll2008h . 

It is clear that the jet could interact with higher density 
gas, and thus be more likely to produce Ni, at early times 
< 1 sec, when the SN shock and MWN are smaller. However, 
the physics at these early times is uncertain: the magnetar 
is likely still contracting to its final radius and the magnetic 
field may even still be growing via a dynamo. Moreover, the 
early time evolution is likely to be more sensitive to the 
details of the explosion mechanism. As described previously 
( fl3.3)L we have verified that the late-time spindown of the 
proto- neutron star and the formation of the collimated jet, 
are not sensitive to the size of the initial cavity we choose 
(which is a proxy for the uncertain SN physics). This will not 
be true of the amount of Ni synthesized. We thus conclude 
that a more careful treatment of the early time contraction 
and spindown of the magnetar, probably together with a 
better understanding of the SN mechanism, are required to 
reliably assess the question of whether the central engine 
responsible for the GRB also contributes significantly to the 
production of Ni during long-duration bursts. 



4 DISCUSSION & IMPLICATIONS 

The calculations presented in this paper, together with our 
previous work (Papers I & II), show that the production of 
a collimated relativist ic jet that can escape the progenitor 
star is a robust consequence of the formation of a B ~ 10 15 
G, P rsj 1 ms neutron star during core-collapse supernovae. 
Paper II shows that such a jet is formed in the low a 
limit in which most of the magnetic energy in the magne- 
tar 's outflow has been converted into bulk kinetic energy - 
a limit motivated by observations of pulsar wind nebulae 
which show that such efficient conversion takes place (e.g., 
iKennel fe Coronit ill 19841 ). In this paper we have focused on 
the opposite limit, that of a highly magnetized outflow in 
which the only conversion of magnetic energy into kinetic 
energy is that which takes place in ideal MHD. Just as in 
the low a limit, we find that a bipolar jet begins to form 
in the first ~ 1 sec after core-collapse and escapes the star 
after ~ 3 — 10 sec (depending on the exact values of B and 
P; see Figs[2][5])- Given the uncertainties in the conversion of 
magnetic energy into kinetic energy in magnetized outflows, 
it is reassuring that the properties of the GRB-producing 
jet are not that sensitive to these details. Physically, this is 
because in all cases the toroidal magnetic field produced by 
the central neutron star builds up in the cavity created by 
the outgoing SN shock, until the field is sufficiently strong to 
drive the flow out in the polar direction; because the speed 
of the SN shock ~ 0.03 — 0.1 c is much less than the speed of 
the magnetar wind ~ c, such a build up of toroidal field is in- 
evitable unless the magnetic energy in the wind is extremely 
small (a < 0.01). 

An important result of our work is that nearly all of the 
spindown energy of the neutron star escapes via the polar 
channel (Fig. [7). There is very little energy transferred to the 
exploding star in the equatorial region. We again find this in 
both the low (Paper II) and high (this paper) a limits. One 
implication of this result is that the central engine powering 
the GRB is unlikely to contribute significantly to energiz- 
ing the SN shock as a whole (although it clearly does so in 



the polar region), at least on timescales > 1 sec after core 
bounce; at earlier times, the dynamics will be sensitive to the 
details of the SN mechanism. We believe that this conclu- 
sion is not specific to the particular central engine considered 
here, but will hold for all magnetically- driven GRB central 
engines. Specifically, we suspect that wide-angle magnetized 
winds from accretion disks, such as are seen in simulation s 
of accretion disks relevant to GRBs (e.g. JProga et al.ll2003h , 
will form a magnetized bubble that will escape via the po- 
lar channel rather than transferring energy to the SN s hock 
as ha s been previ ously hypothesized (Ko hri et al.l l2005) (see 
also (|Aronsll2003h ). 

Observationally, there is a strong association of long- 
duration GRBs with core-collapse SN, in particular en- 
erget i c SN Ic-BL (BL broad line) (IWooslev fe Blooml 
2006: iDella Valid 12006). However, the converse is not true 
(|Soderberg et alJfeOOd : IWooslev fc Blooml l200eSh . Late-time 
radio observations of large samples of Type Ibc SN find 
that even the broad line subset do not show evidence for 
an energetic relativistic outflow, which would produce ra- 
dio emission at late times regardless of whether or not any 
put ative gamma-ray em ission was initially beamed towards 
us jSoderberg et al. 2006). More concretely, Soder berg et alj 
( 2006|) rule out with 84% confidence the hypothesis that ev- 
ery broad- lined SN harbors a GRB. Moreover, the SNe asso- 
ciated with long-duration GRBs are not particularly unusual 
among the class of BL SN in terms of th eir energies , pho- 
tospheric velocities, and Ni masses (e.g., ISoderberg] [2006; 
IWooslev &; Bloomll2006h . A plausible interpretation of these 
data is that some not yet fully understood physics (rota- 
tion?) leads to a class of core-collapse SNe that are un- 
usually energetic and asymmetric (as revealed by spectra- 
polarimetry) , and that produce significant amounts of Ni. 
A subset of those SNe in turn produce relativistic jets that 
power GRBs, but the central engine that powers the GRB 
does not significantly modify the global properties of the 
coincident SNe. Clearly, the detection of a GRB implies 
a preferential viewing geometry, along the rotation axis of 
the system; this is likely accompanied by the detection of a 
modest amount of high velocity eject a and nucleosynthetic 
products f fc|3.4[) . However, there is no indication that the 
global properties or the off-axis appearance of the SNe asso- 
ciated with GRBs are strongly affected by the central engine 
that produces the relativistic jet. Indeed, our results demon- 
strate that after ~ 1 sec, most of the energy produced by a 
spinning-down magnetar escapes along the polar jet rather 
than energizing the SN shock. As argued above, we suspect 
that this result is likely to hold for all magnetically-driven 
GRB models, including disks accreating onto black holes, 
but this remains to be explicitly demonstrated. 

In our calculations we have assumed that ideal MHD 
holds in the magnetized bubble created behind the SN shock. 
It is, however, well known that plasmas with a primarily 
toroidal magn etic field can be unstable to non-axisymmetric 
kink modes (Begclmanl Il998h that would not be captured 
by our axisymmetric simulations (a lthough such instabilitie s 
may be stabilized by rotation; see iTomimatsu et al.1 l200lh . 
This raises obvious questions about the overall stability of 
the magne tic configuration found h ere (s ee also the recent 
results by iMcKinnev fe Blandfordl <|2008h about jet stabil- 
ity). Unfortunately, it is not presently feasible to carry out 
3D simulations of our magnetar wind model with the nec- 
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essary spatial resolution and temporal duration. As noted 
above, however, the production of a relativist ic jet by the 
interaction of a magnetar wind with the surrounding stellar 
envelope appears to be relatively insensitive to the precise 
magnetization of the outflow - as a comparison of the results 
of this paper and those of Paper II demonstrates. The most 
significant difference is in the acceleration of the jet once it 
emerges from the progenitor star: in the highly dissipative 
case (Paper II) the terminal Lorentz factor is larger (~ a 
at the Light Cylinder), while in ideal MHD (this paper), 
the terminal Lorentz factor is somewhat lower, but the out- 
flow remains reasonably strongly magnetized (a > 1) at the 
radii we can study. We thus conclude that unresolved insta- 
bilities are unlikely to change the large-scale evolution of the 
magnetar wind bubble or the resulting production of a col- 
limated jet, but they may well be important in determining 
the acceleration of the jet and the exact ratio of magnetic 
to kinetic energy in the jet at large radii. 

Our results also have implications for the nucleosyn- 
thesis of heavy elements due to the interaction between the 
magnetar- driven jet and the stellar envelope; we largely con- 
firm previous studies on nucleosynthesis in jets and asym- 
metric supernovae, which were based on parameterized mod- 
els for the origin of t he jet / asymmetry (|Maeda et al.l 1 2002 ; 
iNagataki et a l. 2003, 2006). As noted above, most of the en- 
ergy produced by the magnetar is diverted into a jet with an 
opening angle of ~ 5 — 10°. The amount of mass swept-up 
by the jet is ~ O.I M^, most of which is Oxygen in the stel- 
lar model we use (jWooslev et al.ll2002h . Even in our most 
energetic simulations, the temperature at the head of the 
jet is insufficient to produce 56 Ni f 33.41 and Fig. [8}. It is, 
however, sufficient to produce ~ 1O _2 M0 of high velocity 
Ne and Mg, and perhaps a similar amount of high veloc- 
ity Oxygen. Unfortunately, the high velocity Oxygen that is 
observed in some SNe associated with GRBs ( Mazzali et al.l 
2006) can easily be explained as part of the stellar envelope 
that is blown out by the jet, and thus these observations do 
not strongly constrain the properties of the material shock- 
heated by the head of the jet as it propagates through the 
star. The reason that our calculations do not produce a sig- 
nificant amount of 56 Ni is that we initialize our simulations 
roughly 1 second after core bounce, at which time the SN 
shock has moved to ~ 10 9 cm and the post-shock tempera- 
ture at the head of the jet is < 5 x 10 9 K. The temperature 
could be higher at earlier times, but the physics at earlier 
times, and the amount of 56 Ni produced, will likely depend 
sensitively on the physics of the SN explosion itself. 

Our simulations span a factor of ~ 100 in the neutron 
star's spindown power (Fig. [7]). In our high spindown power 
simulations, corresponding to neutron stars with P — 1 ms 
and B ~ 1 — 3 x 10 15 G, we find that the outflow remains 
largely out of causal contact with the surrounding stellar 
envelope. In parti cular, unlike in ID mod els of highly mag- 
netized outflows (|Kennel k Coronit ill 1984 ). in our 2D calcu- 
lations we find that even for high a, the termination shock 
between the wind and the magnetar wind bubble lies outside 
the wind's fast magnetosonic surface, except very near the 
pole (Fig. [I}. Physically, this is because the toroidal mag- 
netic field can escape via the polar jet before it accumulates 
sufficiently in the magnetar wind bubble to drive the termi- 
nation shock to small radii (as occurs in ID). As a result, 
the spindown of the neutron star is essentially identical to 



that of a free wind, i.e., a wind without a surrounding stel- 
lar envelope (Fig. [7|. For our lowest spindown power sim- 
ulation, corresponding to a 3 ms rotator, the termination 
shock does collapse inside th e fast magnetosonic surface (as 
iKomissarov fc Itar kov 2007 found in their 2D simulations). 
Even in this case, however, the neutron star's spindown is 
still quite similar to that produced by a free wind, with at 
most rsj 30% differences (Fig. [7}. In p ractice, these results 
imply that free wind calculations fe.g.. iMetzger et al.ll2007h 
can be usefully used to study the winds from newly formed 
magnet ars. 
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